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ABSTRACT 

It  is  conjectiared  that  for  some  physical  flow  problems  there  may 
be  an  advantage  in  following  energy  exactly  as  one  follows  mass  elements 
in  the  conventional  Lagrangian  formulation  of  fluid  dynamics.  Essential¬ 
ly,  the  argument  is  a  generalization  to  two  dimensions  of  an  idea  due  to 
Enig  who  derived  the  equations  for  one -dimensional  slab  geometry.  The 
equations  are  derived  both  differentially  and  integrally  from  the  Euler i an 
form,  and  they  are  con5)ared  with  the  conventional  Lagrangian  equations. 


NOTATION 

e  specific  internal  energy 

E  total  energy  in  Lagrangian  energy  cell, 

J  Jacobian  of  transformation,  r_z  -  r  z_ 

„  ’  I  Tl  T1  ? 

J-  (r/5)“’J 

P  pressure 

r  radial  coordinate 

t  time 

u  radial  speed 

V  volume 

V  velocity  vector,  (u,0,w) 
w  axial  speed 

z  axial  coordinate 

a  1  or  0,  for  cylindrical  or  Cartesian  plane 
geometries,  respectively 

Tj  Lagrangian  mass  or  energy  coordinate,  z(5,r|,0) 

^  Lagrangian  mass  or  energy  coordinate,  r(5,r],0) 
p  density 

2  2 

iff  total  specific  energy,  ^(u  +w  )+e 
o  subscript  denoting  Initial  values 
A  letter  subscript  denotes  partial  differentiation. 

This  report  derives  and  examines  the  equations 
governing  the  motion  of  energy  elements  in  cylindri¬ 
cal  coordinates  where  symmetry  in  0  is  assumed. 

This  is  analogous  to  the  conventional  Lagrangian 
formulation  of  mass  motion  from  Euler’s  fundamental 
equations.  For  physical  problems  involving  no 


energy  sources,  such  as  impact  problems,  or  for  cer¬ 
tain  problems  wherein  a  time  is  reached  at  which 
there  is  no  further  release  of  energy,  it  is  conjee - 
tiired  that  there  may  be  some  advantage  in  following 
the  flow  of  energy  in  addition  to  mass  flow.  Compu¬ 
tations  of  irr5)act  problems  in  two-dimensional  fluid 
mechanics  in  the  conventional  Lagrangian  formulation 
often  encounter  difficulties  of  mesh  distortion, 
partly  because  of  the  impenetrability  of  mass.  Be¬ 
cause  energy,  however,  can  flow  across  a  collision 
surface,  it  might  be  possible  to  compute  its  motion 
more  smoothly  by  following  such  elements.  Enig^  de¬ 
rived  the  equations  for  one -dimensional  slab  geometry 
and  computed  from  them  the  problem  of  an  infinite 
gas  moving  toward  an  infinitely  rigid  wall.  He 
achieved  good  agreement  with  the  analytic  solution. 

To  my  knowledge,  Enig's  is  the  only  application  of 
this  method. 

We  will  derive  the  equations  differentially, 
leaving  most  of  the  mathematical  details  to  the 
appendixes,  and  will  interpret  the  resulting  equa¬ 
tions  both  differentially  and  integrally. 

We  start  with  Euler's  formulation  of  the  equa¬ 
tions  of  conservation  of  mass,  energy,  and  momentum, 
and  attenpt  to  keep  cylindrical  coordinates  separate 
from  two-dimensional  Cartesian  coordinates 
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0^  +  V.DV  =  0, 


mass. 

(1) 

(oJ). 

T  nvi 

,  +  oJ  — 
t  r 

energy. 

(2) 

r 

+  r 

ril 

[oz^t^o^z^ 

r  momentum, 

(3) 

+  Zp 

or  +0  r 

z  momentum. 

(M 

b 

L  rjt  r|  1 

rjt  n  t 


+  I  OU  +  (du)^  +  (dw)^  =  0, 


(5) 


(d'!()^  +  7  (D'!/+P)u  +  ((Di|r+P)u]^+[(p!|;+P)w]^  =  0, 

(6) 

(ou)^  +  ^  DU^  +  (du^+P)j.  +  (duw)^  =  0,  (7) 

(dw)^  +  ^  pwu  +  (pwu)^  +  (pw^+P)^  =  0.  (8) 

In  the  conventional  way,  we  now  let 

r  -  r(5,Ti,t),  (9) 

z-z(5.Ti,t).  (10) 

Under  this  transformation,  the  partial-derivative 
operators  transform  as 


9r  J  ^  Ti  85  5  8t)/  ’ 

3  i  ('r  .  r  i-') 

j  L  '  8ti  a?y  ’ 


8z 


(11) 

(12) 


where 


J  =  r_z  -  r  z_. 

b  rt  n  ? 


(14) 


Equations  (ll)  through  (13)  are  derived  in  Appendix 

A. 

Applying  Eqs.  (ll)  through  (13)  on  Eq.  {5)>  ve 

have 


p/  +  (r^Z5.-r^z^)o^  +  (r^z^-r^z^)op  +  oj 


+  z  (pu).  -  z.(ou)  +  r  (pw)  -  r  (ow)  =  0. 

'  i  :>  n  q  b  t)  71  'a 


Into  the  above  equation  we  substitute 


P^J  =  (cvT)^  -  p(r^z^^+r^^z^-r^z^^-r^^z^)  , 


and  collect  terms  in  Tj.,  r  ,  z  ,  z..  The  result  is 
b  T]  ri  ? 


This  equation  can  then  be  put  in  the  more  meaningful 
form, 


(pJ)^  +  oJ  -2^  =  rj[p(.z^-w)]^  -  r^[p(z^-w)]^ 

+  z^[R(r^-u)]^  -  z^[p(r^-u)]^.  (15) 


We  note  here  that  Eq.  (15)  is  identically  satisfied 

by 


+  ir(r  z.-r„z  )^  +  (r  z  -r.  z  )^  ,  (13) 
3t  St  t  ^  ^  t^Sr)  U  b  t  ' 


pJ  =  {h°‘  P„(?»Ti). 


rt  =  u, 


z^  =  W, 


(16) 

(17) 

(18) 


because  the  left-hand  side  of  Eq.  (15)  with  the  aid 
of  Eq.  (l6)  is 

(pJ)^  +  oJ  —  =  a(-)  +  an^(^) 

r 


(pj)^  +  dJ  -2^  =  -  p^  a(|)‘"  ^(r^-u) 


Now  looking  back  at  Eqs.  (6),  (7),  and  (8),  we  see 
that  their  forms  are  similar  to  Eq.  (5)^  namely,  a 
term  differentiated  with  respect  to  t,  a  term  in  a 
due  to  radial  motion  in  cylindrical  coordinates,  and 
then  two  terms  differentiated,  respectively,  with  r 
and  z.  We  leave  to  Appendix  B  the  repetitive  alge¬ 
bra  of  applying  Eqs.  (ll)  through  (13)  to  Eqs,  (6) 
through  (8)  to  derive  equations  like  Eq.  (15).  With 
Eqs.  (15)  the  results  are 

(pj)^  +  otJ  rp[o(z^-w)]^  -  r^[p(z^-w))^^ 

+  z  [p(r  -u)]j  -  z  [o(r  -u)]  mass,  (I9) 

T]  ^  T)  b  ^  n 

(pjilf)^  +  aj(pt!t+p)^=  r^[pi|r(z^-w)-Pw]^-r^[plt(z^-w)-Pw]^ 

+z  [pt(r  -u)-Pu]  -z  [D.l.(r  -u)-Pu]  energy,  (20) 

T]  U  S  b  ^  q 
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2 

(oJu)^  +  aJ  -2^=  r^(pu(z^-w)]^-r^[Du(z^-w)]^ 

+  2  [pu(r  -u)-p]  -z  [pu(r,-u)-P]  r  momentum. 
T\  X  .  b  9  ^  Ti 

(21) 

and 


ventional  Lagrangian  formulas  in  two  dimensions. 

To  follow  energy  elements,  instead  of  mass 
elements,  let  us  ask  what  are  the  consequences  of 
satisfying  Eq.  (20)  identically.  To  accoir5)lish  this, 
let 


(r>Jw)^  +  aJ  2^=  r^[pw(z^-w)-p]^-r^[pw(z^-w)-p]^ 


pJ'lf  =  pQ(?,r|)i|t^(?,Ti)(|)“'  . 


(27) 


-  z  [pw(r  -u)]p-2p.[pw(r  -u)]  z  momentum.  (22) 

T|  ^  -3  ■=>  XT) 

Let  us  consider  what  has  been  done.  We  started  with 
the  four  Eulerian  equations,  (5)  through  (8),  in  the 
five  dependent  variables,  o^u,  w,  P,  and  e,  and  the 
three  independent  variables,  r,  z,  and  t.  We  assume 
that  we  also  have  an  equation  of  state  relating  p, 

P,  and  e  so  that  we  have  five  equations  in  five  un¬ 
knowns.  With  the  transformations,  Eqs.  (9)  and  (lO), 
We  have  created  two  more  dependent  variables,  r  and 
z;  hence,  Eqs.  (19)  through  (22),  complemented  with 
an  equation  of  state,  constitute  five  equations  in 
seven  unknowns  in  the  new  independent  variables 
X],  and  t.  In  the  conventional  Lagrangian  fornnila- 
tion,  Eq.  (I9)  is  identically  satisfied  by  replac¬ 
ing  it  with  Eqs.  (16),  (17),  and  (I8).  Thus,  we 
have  seven  equations  in  the  seven  unknowns,  r,  z. 


u,  w,  0,  P, 

and  e.  With  the  aid 

of  Eqs.  (16)  and 

(17),  we  note  that 

(oJ0)t 

+  aJD0  ^  = 

■  <l>Vt 

-  “•(|)°'Po  7(^t" 

d)^ 

or 

(oJ0)^ 

+  aJD0  ^  = 

> 

(23) 

where  0  is  any  dependent  variable.  Thus  with  ] 

Eqs. 

(23),  (17), 

and  (18), 

Eqs.  (20), 

(21),  and  (22) 

sin5)lify  to: 

0  — 1 

oL  t  orj 

=  (1)  l-r 

energy. 

(24) 

f  r \0c  r 

P-+Z,  P  ] 

?  1  T1 

r  momentum, 

(25) 

and 

II 

0 

c 

(T)‘"[-r, 

b 

P  +r  P,] 
n  n  b 

z  momentum. 

(26) 

The  above  three  equations  are  recognized  as  the  con- 


Then, 


(pJilt)t  =  -Pj,  iIIq  »(!)“■ 


(28) 


and  the  left-hand  side  of  Eq.  (20)  becomes 

-P  !|r  ct(-^)®  ifr  -u- 2^]  .  (29) 

^o  *^0  'r''  rL  t  p\ifj  ' 

Examining  the  right-hand  side  of  Eq.  (20),  we  see 
that  it  also  will  equal  zero  if 


(30) 


(51) 


Hence  energy  flows  with  speeds  (1+  P/p(jf)u  and 
(1+  p/pi!r)w  in  the  radial  and  axial  directions,  re¬ 
spectively.  It  should  be  emphasized  that  all  our 
dependent  variables  except  r  and  z  retain  their  same 
physical  meaning  as  in  the  conventional  Lagrangian 
formulation;  specifically,  u  and  w  are  still  the 
mass  particle  speeds.  However,  through  their  de¬ 
pendence  on  the  new  independent  variables  9  and  n, 
they  are  correlated  with  the  initial  position  of  a 
given  energy  element.  For  impact  problems  in  which 
one  can  assume  that  initially  there  is  only  kinetic 
energy  of  some  projectile  or  moving  medium,  the 
initial  position  of  a  given  energy  element  will  be 
the  same  as  the  initial  position  of  the  correspond¬ 
ing  mass  element. 

To  complete  the  formulation  of  our  problem,  we 
must  now  ask  what  are  the  consequences  of  applying 
Eqs.  (27),  (30),  and  (31),  to  Eqs.  (I9),  (21),  and 
(22).  Because 


(pj)^  =  ir^[(pJ']r)^-pJ'|rt]  , 


and  with  Eq.  (28)  the  left-hand  side  of  Eq.  (I9) 
becomes 


3 


Then  substituting  Eqs .  (50)  and  (51)  into  the  above 
and  the  right-hand  side  of  Eq.  (I9),  we  have 


where  V  is  the  moving  volume  element  enclosed  by  the 
moving  surface,  S.  By  convention,  dS  points  outward 
from  the  given  energy  element.  The  right-hand  side 
is  then  the  mass  flux  into  a  particular  energy  element; 
thus,  the  operator,  d/dt,  on  the  left-hand  side  im¬ 
plies  following  the  given  energy  element  in  the  fixed 
frame.  In  other  words,  Eqs.  (50)  and  (51)  imply  that 
the  rate  of  growth  of  the  volume  of  a  given  energy 
element,  i.e.. 


roi  ^  ' 

-r  (Pwif  ^).  +  z  (Pu'ir  ^)_-  z^(Pu!!/  ^)  mass  . 


(32) 


is  such  that 


In  Appendix  C  we  repeat  this  procedure  for  Eqs,  (21) 
and  (22).  The  results  are 


(36) 


■^]  "  rj,[Pwu,r^]  -r  [Pvm.r^] 

TPi’lf  ^  'll 


+  z  [pu^-lf  ^-P  ] --2^  [Pub'll  "^-p]  ,  r  momentum,  (55) 


We  think  of  the  left-hand  side  of  Eq.  (56)  as  ex¬ 
pressed  in  Lagrangian  energy  coordinates,  (5,r|,t), 
and  of  the  right-hand  side  as  a  function  of  (r,z,t). 
Further,  the  operator,  d/dt,  in  Eq.  (55)  is 


+z  [puwijr  ^l^-z.[Puwi!f  ,  z  momentum.  (3^) 

“n  ^  "j  “n 


=  (tt-)  +  (1+  “)v-V 

dt  ^Bt^r,z  ^  p'y 


(37) 


and  the  transformation  of  this  operator  to  Lagrangian 
energy  coordinates  is  simply 


Because  Eq.  (52)  originally  came  from  the  mass 
conservation  equation  in  Eulerian  form,  it  should 
express  mass  conservation  as  one  follows  in  time  a 
particular  energy  element  or  cell,  which  is  identi¬ 
fied  by  its  original  position  through  (§,1^).  To 
show  that  this  is  indeed  the  case,  we  will  rederive 
Eq.  (52)  integrally  with  the  hope  that  its  physical 
content  will  be  more  plausible.  Equations  (50)  and 
(51)  tell  us  that  energy  moves  with  radial  and  axial 
speeds,  respectively,  as 


dt 


(— ) 


(38) 


Equations  (56)  and  (57)  are  derived  in  Appendix  D. 

In  the  conventional  Lagrangian  coordinates,  wherein 
mass  is  followed,  the  equations  analogous  to  Eqs. 

(56)  and  (57)  are 

=  V.^  ,  (39) 

J' 

and 


(1+  “)u  and  (1+  — )w  . 
pir  P’lf 


d  .a  V 

dt  ^Bt^r,z 


(^0) 


Hence,  the  velocity  of  a  mass  particle  relative  to 
the  moving  energy  element  is 


7  -  (1.  ^)7  = 


Therefore,  to  express  mass  conservation,  we  write 


To  return  to  Eq.  (55),  we  first  rewrite  is  as 


We  then  transform  Eq.  (4l)  to  Lagrangian  energy 
coordinates , 


pdV 


P  -♦ 

-  7  V  dS  , 

ij/ 


(35) 


(— ) 


(v.|7)j'dVo  , 


{h2) 
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because  dv(r,z)  =  J'dV^(5,Ti)  =  J'd^drj.  The  bar  over 
V  reminds  us  that  this  operator  is  to  be  trans¬ 
formed  to  (^,r|)  coordinates.  We  now  can  take  the 
operator,  (^/dt)^  ^  on  the  left-hand  side  under  the 
integral, 

r(c^')tdVo=  r  (v-f  v)j-dv^ ,  (43) 

or 

r  =  I  [(pu,r^)^  +  (pwr^)^ 

+  |(Pu'r^)]j'dV^  .  (44) 


/  ^  \  d 

"  dt  ' 

Because  the  second  and  third  terms  can  be  combined 
as 

v*d(1+  ^)v  =  V.pv  +  v4  V  ,  (49) 

pi(f  ^!/ 

we  see  that  we  are  left  with  Eq.  (l). 

The  above  arguments  of  plausibility  are  easily- 
extended  to  the  two  momentum  equations,  Eqs.  (33) 
and  (3^).  For  the  r  momentum,  the  equation  of  con¬ 
servation  is 


In  the  integrand  on  the  left-hand  side  of  Eq. 
we  substitute  from  Eq.  (27)^ 

pJ^  =  0^  i.!fo  ^  ^  (27a) 

and  carry  out  the  transformations  indicated  on  the 
right-hand  side  with  Eqs.  (11)  and  (12), 

I  =  f  { 

O  O 

+  r^(PwJ.-^)^  -  r^(Pw,r^)  J(|)°'+  ^  Puili'^  ^']^o  ' 

Because  our  volume  element,  dV^,  is  arbitrary,  the 
integrands  in  Eq.  (^5)  tnay  be  equated,  and  we  recog¬ 
nize  that  we  have  Eq.  (32),  after  elimination  of  J‘ 
in  the  last  term  with  Eq.  (27a). 

Finally,  we  can  make  one  other  check;  we  can 
go  from  Eq.  (^3)  back  to  the  fixed  Eulerian  frame 
with  the  help  of  Eqs.  (36)  through  (38)  to  see  wheth¬ 
er  we  get  Eq.  (l).  We  first  differentiate  the  inte¬ 
grand  on  the  left-hand  side  of  Eq.  (^3)^ 

I  (o^J'+^’^)dV^  =  I  (V-£  v)j'dV^  ,  (46) 

''o  ^o 

and  eliminate  J‘  using  Eq.  (36), 

{  =  I  (^^7  7)j'dv^  .  (47) 

o  o 

We  now  transform  to  the  fixed  Eulerian  frame,  remem¬ 
bering  that  for  p^. 


r  Dudv  =  .^  pu[v-(l+  f^)7]-S  -  ^  P(dS)'',  (50) 

V(t) 

where  the  first  integral  on  the  right-hand  side  is 
the  r  momentum  flux  into  the  energy  element,  and  the 
second  integral  is  the  net  force  in  the  r-direction. 
Transforming  to  volume  integrals  on  the  right-hand 
side,  we  have 


""  l(t) 


pudV  = 


~  v  dV 


dV  . 


(51) 


Again,  transforming  to  Lagrangian  energy  coordinates, 
we  have 


I  (ouj’)tdv„  =  I 

^o  o 

+  I  Pu2,|,-1  .  p^ll’dV^  .  (52) 


Eliminating  J'  through  Eq.  (27a)  in  the  integrand 
on  the  left-hand  side  and  in  the  third  term  on  the 
right-hand  side  and  carrying  out  the  indicated  trans¬ 
formation  using  Eqs.  (ll)  and  (12)  on  the  right-hand 
side,  we  have 

I  =  I  {  [^^(PuS-h^- 

o  ^o 

+  r^(Puwii'^)^  -  r^(Puwir^) ^(1)“- 

+  —  -  [z  P.-ZpP  ](5)°'|dV^  . 

ro  ^  T]  s  ?  ^  ?  >1  o 

(53) 
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Because  dV^  is  arbitrary,  Eq.  (53)  Is  Eq.  (33) » 

To  see  if  Eq.  (52)  yields  Eq.  (3)  when  trans¬ 
formed  back  to  a  fixed  Eulerian  frame,  we  first  dif¬ 
ferentiate  the  integrand  on  the  left-hand  side, 

I  [(Du)^J'+puJ'^]dV^  =  I  vJdV°-|  JdV°  .  (5'^) 

^  o  ^  o  o 

Substituting  Eqs.  (36)  through  (38)  and  transforming 
to  (r,z,t)  coordinates,  we  have 

=  r  7  dV  -  r  P  dV  .  (55) 


The  second  and  third  terms  in  the  integrand  on  the 
left-hand  side  may  be  combined  as 


v*ou(l+  ^)v  =  V»ouv  +  V 

^  O’J/ '  )(/ 


Hence,  Eq.  (55)  becomes 

r  [(nu)  +  V*Duv]dV  =-[  P  dV  ,  (56) 

V(t)  ^  V(t)  ^ 

which  is  Eq.  (3)  for  arbitrary  dV.  The  argument 
for  Eq.  (3^),  the  z-momentum  equation,  follows  iden¬ 
tically  to  the  above  argument  for  the  r  momentum. 

Thus,  Eqs.  (32)  through  (3^)  are  the  fundamental 
differential  equations  to  be  solved  in  order  to  fol¬ 
low  energy  elements.  The  terms  should  be  elim¬ 
inated  in  finite  differencing  because  there  will  be 
cells  with  no  total  energy  in  them.  Some  siiE5)lifi- 
cation  can  be  achieved  by  subtracting  from  each  mo¬ 
mentum  equation  the  corresponding  speed  times  Eq.(32). 
Taking  Eq.  (33)  and  subtracting  from  it  u  times  Eq. 

( 32) ,  we  have 


D  ’if  (*^)°''  =  r-u  Pwt|("^  -  r  u_Pw^  ^ 

^o  o^r'  t^  §  rt  ^  11  5  ^ 


+  z^(u^r^-P5>z^u^Pur^-P^), 


D  t  u .  =  (^)nPw(r.u  -r  u-)+Pu(z  u.-z  u  ) 
■o-'o  t  L  ^  ^  T1  ri  T)  §  P  ri'' 


[z  P,-Z,P  )  . 

'  Ti  5  ?  n'J 


Similarly,  subtracting  w  times  Eq.  (32)  from  Eq.  (3^), 
we  have 

PoVt  ■ 

Finally,  we  can  eliminate  'If  ^  in  Eq.  (32)  by  carrying 

out  the  differentiations  and  multiplying  through  by 
2 

ilf  .  The  resulting  equation  is. 


Equations  (57)  and  (58)  should  be  compared  with 
Eqs.  (25)  and  (26),  their  analogues  in  the  conven¬ 
tional  Lagrangian  formulation;  similarly,  Eq.  (59) 
should  be  con5)ared  with  Eq.  (2^).  Certainly  the 
equations  in  Lagrangian  energy  coordinates  are  more 
con5>licated  and  their  computations  will  require  more 
machine  time.  When  differenced,  each  acceleration 
equation  will  have  three  times  as  many  factors  as 
its  conventional  Lagrangian  analogue.  Because  we 
are  following  energy  cells,  it  would  seem  reasonable 
when  differencing  these  equations  to  consider  as 
a  cell-centered  variable;  hence,  we  ascribe  a  value 
for  u  and  w  at  the  energy  cell  center  by  averaging 
the  values  at  the  four  corners.  Further,  because 
all  the  equations  involve  velocity  gradients,  it 
would  make  additional  sense  to  have  cell-centered 
velocities.  If  one  examines  some  of  the  more  popu¬ 
lar  methods  of  finite  differencing  the  conventional 
formulas  and  extends  these  methods  to  Eqs.  (57) 
through  (59) ^  he  will  note  that  the  mass  accelerations, 
u.  ,  w  ,  of  a  vertex  point  in  the  energy  mesh  depend 
through  (|f  on  information  within  the  four  surrounding 
Lagrangian  energy  cells  as  well  as  at  each  corner  of 
these  cells.  To  illustrate,  we  give  the  explicit 
contribution  to  the  r  acceleration  of  cells  whose 
midpoints  are  1  and  4  in  the  following  figure. 

That  is,  we  difference  Eq.  (57)  by  assuming  the  ac¬ 
celeration  in  the  r-direction  of  point  0  to  be  made 
up  of  four  parts  contributed  by  pairs  of  cells,  1 
and  4,  2  and  3>  1  and  2,  and  4  and  3-  The  contri¬ 
butions  by  cells  1  and  k-  and  2  and  3  are  averaged, 
as  are  the  contributions  by  cells  1  and  2  and  4  and 
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3  (see  Ref.  2,  p.  l8) .  The  result  is 

Cr^f  (^5)“ 

*  E2+E3  •  •  •  "^  Ej^+E2  •  •  ■  ■^  Ej+Ej^  •  •  •  > 

where  ?g  =  Ur^+rJ  or  Kr^+r^^)  and  E^^  = 

where  V  is  the  volume  at  t  =  0  for  cell  i.  In  con- 
o 

trast,  the  same  process  for  Eq.  (25)  yields 

plus  three  similar  terms,  where  m^  = 
have  assumed  that  all  dependent  variables  are  known 
in  some  average  way  at  all  the  lettered  and  numbered 
points  in  the  figure.  In  practice,  this  could  be 
achieved  for  points  a,  b,  c,  d  by  getting  information 
from  cells  surrounding  those  depicted,  or  such  terms 
could  be  ignored. 

Interfaces  between  different  materials  will  pose 
an  additional  complication  here,  because  we  will  have 
to  know  them  in  order  to  choose  the  appropriate  equa¬ 
tion  of  state;  i.e.,  in  an  energy  cell  that  overlaps 


two  or  more  mass  cells,  there  will  be  contributions 
to  the  cell-centered  pressure,  density,  and  specific 
total  energy  from  each  of  the  mass  constituents. 

Mass  interfaces  can  be  found,  however,  because  u^, 
w^  ere  the  mass  accelerations  of  a  vertex  point  in 
the  energy  mesh.  Ve  would  have  to  solve  two  addi¬ 
tional  equations. 


^tt  =  ^  ^tt  =  ''t 


for  instance,  to  follow  mass  points  and  then  have  to 
test  r  and  z  against  these  values  to  find  the  loca¬ 
tion  of  the  energy  points  relative  to  the  mass  points. 
Certainly,  it  would  be  prudent  to  first  try  this 
method  on  a  single -material  problem. 

Some  remarks  on  the  order  of  solving  these 
equations  are  warranted,  because  this  order  demon¬ 
strates  an  additional  difficulty  not  encountered  in 
the  conventional  Lagrangian  formulation.  A  proposed 
order  of  solution,  assuming  that  all  dependent  vari¬ 
ables  except  u  and  w  are  known  at  n,  and  that  u  and 
w  are  known  at  n  -  J,  is: 

1.  Solve  Eq.  (59)  ^ 

n^.1 

2.  Solve  Eq.  (57)  for  u  ^  , 

3.  Solve  Eq.  (58)  for  . 


From  these  three  equations,  we  can  get 


1/ )  i^\  /  2  2Nn+| 

e  ^  =  ^(i|;  +t  )  -  (u  ^  . 


5 .  Solve  0^ 


P  t 
o  ’  o 


(J’) 


n  n+1 


6.  Assuming  a  Mie-Gruneisen  equation  of  state, 
solve  for. 


P*  .  f^(o*)  +  . 


7. 

p^ilr 

8.  z*-z"  =  (At)"'^^(^l+  ’ 

P*i(r  ^ 

where 

(At)"^  =  -  t"  . 


From  r*,  z*  we  now  find  a  new  value  for  J'  to 
substitute  into  Eq.  (5)  and  iterate  on  Eqs.  (5) 
through  (8)  until  consistency  is  attained  between 
J',  r^,  z*.  Call  the  final  values  r^'*’^, 


T 


n+i  n+i 

0  ,  n  .  Finally  solve  for 


and  the  inverse. 


=  p"  +  ?P"^,  =  o"  +  2o“^^  . 

Thus,  all  quantities  are  advanced  to  n  +  1  except  u 

and  V  which  are  advanced  to  n  +  ^. 

Comment  from  any  reader  who  has  used  this  method 

in  a  two-dimensional  code  is  invited. 
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APPENDIX  A.  DERIVATION  OF  EQS.  (ll),  (12),  and  (13) 
Consider  the  transformation  from  (r,2,t)  to 
and  the  inverse, 

^  ^  =  5(r,z,t) 

z  =  z(5>Ti/ti’)  n  =  'n(r,z,t) 
t  =  t(0,0,t’)  t’  =  t’(0,0,t)  , 

such  that  t  =  t'.  We  retain  the  prime  on  t  associated 
with  the  n,  t'  set  of  independent  variables  for 
convenience.  The  operators  of  partial  differentiation 
are  related  by 


dr  \  an  a§  ’ 


2-  =  „  9.  +  .  Sl. 
dz  ^2  an  ^2  a§  ’ 

^  =9-  +  ^  ^  +  F  9. 

3t  Bt’  t  ^rj 


We  think  of  the  coefficients  in  the  above  equations, 
r\^,  and  5^,  as  functions  of  r,  z,  t, 
so  we  must  find  the  transform  of  these  functions  in 
order  to  get  the  operators  on  the  right-hand  side  as 
functions  of  our  new  independent  variables  n,  t. 

To  accomplish  this,  we  consider  how  the  differentials 
transform. 


dri  =  n 


Now  we  can  invert  Eq.  (A-2)  symbolically. 


<3-n  h  (cofactor  of  a^  .  transposed)  dz  ,(A-4) 


where  a. ,  is  the  matrix  in  Eq.  (A-2)  and  la.  .|  is 
the  determinant  of  a^^.  Explicitly,  Eq.  (A-4)  is 


z  -r  r  z, -r,  z  \  /  dr 
U  Ti  11  t  t  n  \  / 


-  |a.  .(^  .  -z^  .  r,  -(r^Z^-r^z.)  dz  ..  {A-5) 


0  0  r-z  -r  z  /  \ dt 

^  5  ri  ri  ?/  \  . 


Equations  (A-3)  and.  (A-5)  demand  that  the  correspond¬ 
ing  elements  in  their  respective  matrices  be  equal. 

In  addition,  we  see  that 


=  J 


our  previous  notation.  Hence, 


Tij,  =  -J  ^2^,  =  J  = '-J'^(r^z^-r^z^). 


Substituting  Eqs.  (A-6)  into  Eqs.  (A-l)  we  have  Eqs. 

(11),  (12),  and  (15)  in  the  body  of  this  report. 

APPENDIX  B.  DERIVATION  OF  EQS.  (20,  (21)  and  (22) 
Applying  Eqs.  (11),  (.12),  and  (13)  toEq.  (6) 
yields 

( ot)yt(r.j.z^-r^z;|.)(  pil,)^+(r^z.|.-r.|.z^)(  oi)^+aJ(nt+P)^ 

+2  [(Di|r+P)u]  -z^[(nt+P)u]  +r,[(o'!-+P)w]  -r  [( ot+P)v] .  =  0. 

(B-1) 

Into  Eq.  (B-1),  substitute 

(P>^)y  =  ,  (b-2) 


8 


The  result 


and  then  collect  terms  in  r_,  r  ,  z  ,  z.. 

?  n’  Ti'  ? 


-  r^{p'J'Z^^+(p\lf)^z^-((Dtl(+P)w)^j- 
+  z^|o<»r^^+(oilf)^r^-[(pil(+P)u]^| 


APPENDIX  C.  DERIVATION  OF  EQS.  (33)  AND  (3^) 

With  Eq.  (27),  the  left-hand  side  of  Eq.  (21) 
becomes 

!2^  (S)<^ .  +  r_2;^  (S)<^i  u  +  a  -°^S)“  3i- ,  (c-i) 

i(r  r  t  Lit  r  J.  il[^r  r^'  ' 


which  is  Eq.  (20). 

Repeating  this  procedure  on  Eq.  (7), 

2 

( ou)^J+(r^z^-r^z^)(  pu)^+(r^z^-r^z^)(  Du)^4aJ 

+z  (ou^+P).-Zp(pu^+P)  +r,(Duw)  -r  (puw)p  =  0;  (B-4) 

^  S-D  'IS  "a 


Then,  with  Eqs.  (30)  and  (51),  Eq.  (21)  is 

2 

P  ^  (*^)^  ‘If  Tu. -ilf  ^^.u  -  aiff  ^  £H_~| 
^o^o^r'  L  t  t  roJ 


(pu)^J  =  (puJ)^  -  ou  ; 


2 

(nuJ)^+aJ  r^[puz^^+(nu)^z^-(ouw)J 

-(pu^+P)  J 

which  is  Eq.  (21) . 

Finally,  repeating  this  procedure  on  Eq.  (8), 

( ow)^J+{r^z^-r^z^)(  „w)^+(r^z^-r^z^)(  pw)^+aJ 

+Z  (pwu).-Zp(owu)  +r.(pw^+p)  -r  (pw^+P),  =  0  j  (B-7) 

^  S  -D  >111  s 

(dw)^J  =  (pvj)^  -  pwj^  ;  (B-8) 

(pwj)^4txj  r^[pwz^^+(Dw)^z^-(pw2+P)J 

P^^+P)  P'')  ft 

-  (dwu)|] 

which  is  Eq.  (22). 


=:  r^(Pwuir^)^  -  r^(Pwu?r^)p 

+  z  (pu^.r^-p)^  -  z^(pu^,r^-p) 

“H  Di  Tn 


which  is  Eq.  (35). 

Repeating  the  above  procedure  on  Eq.  (22),  the 
left-hand  side  is 


-  Po^[f(|)"-?^-|(|)"> 


P  ’!/ 

'^O  0/  wu 


D  'J;  ff  , 

O  O/gNtzl 


-  -T • 


Then,  with  Eqs.  (50)  and  (51),  Eq.  (22)  is 


-r  (Pw^t!r**^-P)p  +  z  (Puwlf‘^)  -  z  (Puwiir”^)  ,  (C-7) 

'I  S  T]  b  ^5  T\ 

Which  is  Eq.  (5^) . 

APPENDIX  D.  DERIVATION  OF  EQS.  (36)  AND  (37)  FROM 
EQS.  (30)  AND  (31) 

By  definition. 
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r 

(D-l) 

+ 

'3  Ti 

r,  r 
s  ^ 

+  (^w)^ 

•?  n 

Zp  Z 
b  n 

J’  =  ijf 


r.  r 
=S  ri 
z„  z 

V 


(D-2) 


Differentiating  Eq.  (D-2)  with  respect  to  t, 


j;  =  (f)" 


r-*  r  4. 
■%t  nt 


-  a(?r^  ^ 


(D-3) 


Into  Eq.  (D-3)  we  now  substitute  Eqs.  (30)  and  (51), 


j;  =  t(0u)  +  ^  0u  +  (0w)  Jj'  , 


+  a^J-  , 

r  ' 

(D-10) 

(D-11) 


which  is  Eq.  (36). 

For  Eq.  (37),  we  first  rewrite  Eq.  (13)  in  the 
body  of  report  in  the  form,  using  Eqs.  (11)  and  (12), 

(L.)  =  (^)  _  (^)  (^)  -  (15)  (1_) 

'■9t^r,z  '■at^5,Ti  'srbjt  '3t^§,r| 

(D-12) 


r^=(l+^)u,  z^=(l+^)w  , 


and  for  convenience  let  l  +  —  =  0.  Then, 

P’ii  ^ 


Jt  = 


(0u)  (j^u) 


+  (f)'" 


(0w)r  (1^) 


(g);„  ■ 

(D-13) 

Then,  substituting  Eqs.  (30)  and  (31)  into  the  above, 
we  have 


^  /Tta-l  0u 
+  a(^)  ^  J  . 


(D-4) 


(■r:r)ff  “  >  (D-14) 


Now  consider  0u  and  0w  as  functions  of  r  and  z;  then,  which  is  Eq.  (37). 


le  chain  rule, 

(flu).  =  r^(0u)  +  Z  (0u)  ; 

(D-5) 

(isu)  =  r  (^u)  +  z  (0u)  ; 

r)  y\  i.  yy  ^ 

(d-6) 

(0w),  =  r,(0w)^  +  z^(0w)2  J 

(D-7) 

(0w)^  =  r^(0w)^  +  z^(S^w)^  . 

(d-8) 

Substituting  the  above  relations  into  Eq.  (D-4),  and 
in  the  last  term  substituting 


J  =  J7(f)^  , 


(D-9) 


we  have 


j;  =  (j)""  (0u) 


Zp  Z 
S  r 


(0^), 


Zp  z 

?  n 

Z-  Z 
?  T1 
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